{ float angles[90] = { 0.500000 ,1.500000 ,2.500000 ,3.500000 ,4.500000 ,5.500000 ,6.500000 ,7.500000 ,8.500000 ,9.500000 ,10.500000 ,11.500000 ,12.500000 ,13.500000 ,14.500000 ,15.500000 ,16.500000 ,17.500000 ,18.500000 ,19.500000 ,20.500000 ,21.500000 ,22.500000 ,23.500000 ,24.500000 ,25.500000 ,26.500000 ,27.500000 ,28.500000 ,29.500000 ,30.500000 ,31.500000 ,32.500000 ,33.500000 ,34.500000 ,35.500000 ,36.500000 ,37.500000 ,38.500000 ,39.500000 ,40.500000 ,41.500000 ,42.500000 ,43.500000 ,44.500000 ,45.500000 ,46.500000 ,47.500000 ,48.500000 ,49.500000 ,50.500000 ,51.500000 ,52.500000 ,53.500000 ,54.500000 ,55.500000 ,56.500000 ,57.500000 ,58.500000 ,59.500000 ,60.500000 ,61.500000 ,62.500000 ,63.500000 ,64.500000 ,65.500000 ,66.500000 ,67.500000 ,68.500000 ,69.500000 ,70.500000 ,71.500000 ,72.500000 ,73.500000 ,74.500000 ,75.500000 ,76.500000 ,77.500000 ,78.500000 ,79.500000 ,80.500000 ,81.500000 ,82.500000 ,83.500000 ,84.500000 ,85.500000 ,86.500000 ,87.500000 ,88.500000 ,89.500000 }; float all[90] = { 0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.157688 ,0.480366 ,0.271876 ,0.202835 ,0.388946 ,0.344739 ,0.305480 ,0.300315 ,0.292566 ,0.288863 ,0.284759 ,0.288427 ,0.272461 ,0.256115 ,0.253577 ,0.252343 ,0.251450 ,0.251050 ,0.256947 ,0.260073 ,0.274027 ,0.277225 ,0.284345 ,0.283710 ,0.283365 ,0.283367 ,0.286163 ,0.288472 ,0.288751 ,0.289256 ,0.290034 ,0.291031 ,0.288455 ,0.269891 ,0.199981 ,0.130150 ,0.127476 ,0.124950 ,0.122544 ,0.120275 ,0.118122 ,0.116089 ,0.114136 ,0.112292 ,0.110529 ,0.108873 ,0.107295 ,0.105787 ,0.104358 ,0.102985 ,0.101684 ,0.100442 ,0.099276 ,0.098167 ,0.097105 ,0.096074 ,0.095112 ,0.094205 ,0.093343 ,0.092513 ,0.091723 ,0.090978 ,0.090276 ,0.089609 ,0.088972 ,0.088379 ,0.087822 ,0.087297 ,0.086801 ,0.086339 ,0.085912 ,0.085511 ,0.085135 ,0.084798 ,0.084477 ,0.084196 ,0.083929 ,0.083712 ,0.083497 ,0.083322 ,0.083170 ,0.083040 ,0.082939 ,0.082867 ,0.082811 ,0.200259 }; float vxd[90] = { 0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.005548 ,0.016283 ,0.017250 ,0.018439 ,0.023814 ,0.024806 ,0.024193 ,0.030709 ,0.028187 ,0.023710 ,0.022023 ,0.021365 ,0.020795 ,0.020492 ,0.026275 ,0.029194 ,0.028544 ,0.026466 ,0.024867 ,0.024310 ,0.023808 ,0.023432 ,0.025625 ,0.027122 ,0.026393 ,0.025677 ,0.025030 ,0.024407 ,0.023836 ,0.023279 ,0.022759 ,0.022280 ,0.021822 ,0.021397 ,0.020986 ,0.020597 ,0.020225 ,0.019877 ,0.019540 ,0.019233 ,0.018926 ,0.018633 ,0.018374 ,0.018113 ,0.017868 ,0.017628 ,0.017406 ,0.017185 ,0.016996 ,0.016811 ,0.016629 ,0.016442 ,0.016280 ,0.016129 ,0.015984 ,0.015843 ,0.015705 ,0.015576 ,0.015458 ,0.015342 ,0.015229 ,0.015127 ,0.015032 ,0.014945 ,0.014859 ,0.014778 ,0.014709 ,0.014640 ,0.014572 ,0.014519 ,0.014459 ,0.014414 ,0.014363 ,0.014336 ,0.014291 ,0.014264 ,0.014240 ,0.014215 ,0.014198 ,0.014188 ,0.014174 ,0.014171 }; float ftd_sit[90] = { 0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.157688 ,0.480366 ,0.271876 ,0.201369 ,0.183432 ,0.110834 ,0.060687 ,0.060247 ,0.056245 ,0.055615 ,0.053971 ,0.059605 ,0.045192 ,0.030049 ,0.028398 ,0.027778 ,0.027248 ,0.026989 ,0.032818 ,0.035683 ,0.049195 ,0.051774 ,0.058110 ,0.056531 ,0.055079 ,0.053833 ,0.055207 ,0.055947 ,0.054521 ,0.053149 ,0.051888 ,0.050696 ,0.048750 ,0.045994 ,0.044971 ,0.044018 ,0.043114 ,0.042265 ,0.041453 ,0.040684 ,0.039954 ,0.039266 ,0.038603 ,0.037987 ,0.037387 ,0.036818 ,0.036294 ,0.035781 ,0.035298 ,0.034830 ,0.034390 ,0.033963 ,0.033578 ,0.033206 ,0.032846 ,0.032490 ,0.032166 ,0.031863 ,0.031573 ,0.031294 ,0.031025 ,0.030771 ,0.030535 ,0.030309 ,0.030090 ,0.029889 ,0.029701 ,0.029526 ,0.029357 ,0.029199 ,0.029058 ,0.028922 ,0.028792 ,0.028682 ,0.028569 ,0.028477 ,0.028382 ,0.028317 ,0.028238 ,0.028181 ,0.028131 ,0.028085 ,0.028051 ,0.028029 ,0.028006 ,0.027998 }; float tpc[90] = { 0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.157688 ,0.480366 ,0.271876 ,0.202835 ,0.388946 ,0.336203 ,0.263293 ,0.257983 ,0.250078 ,0.246202 ,0.241913 ,0.245379 ,0.229199 ,0.212620 ,0.209835 ,0.208338 ,0.207165 ,0.206467 ,0.212047 ,0.214840 ,0.228443 ,0.231266 ,0.237993 ,0.236944 ,0.236162 ,0.235709 ,0.238019 ,0.239819 ,0.239569 ,0.239513 ,0.239699 ,0.240079 ,0.242858 ,0.269352 ,0.177773 ,0.108412 ,0.106184 ,0.104082 ,0.102078 ,0.100187 ,0.098393 ,0.096700 ,0.095072 ,0.093539 ,0.092069 ,0.090688 ,0.089375 ,0.088118 ,0.086928 ,0.085784 ,0.084700 ,0.083663 ,0.082695 ,0.081772 ,0.080887 ,0.080027 ,0.079226 ,0.078471 ,0.077753 ,0.077063 ,0.076403 ,0.075783 ,0.075198 ,0.074643 ,0.074111 ,0.073617 ,0.073154 ,0.072716 ,0.072303 ,0.071918 ,0.071563 ,0.071229 ,0.070915 ,0.070635 ,0.070367 ,0.070133 ,0.069910 ,0.069731 ,0.069551 ,0.069405 ,0.069279 ,0.069170 ,0.069086 ,0.069027 ,0.068979 ,0.186431 }; float set_etd[90] = { 0.000000 ,0.000000 ,0.000000 ,0.000000 ,0.157688 ,0.480366 ,0.271876 ,0.202835 ,0.388946 ,0.344739 ,0.305480 ,0.300315 ,0.292566 ,0.288863 ,0.284759 ,0.288427 ,0.272461 ,0.256115 ,0.253577 ,0.252343 ,0.251450 ,0.251050 ,0.256947 ,0.260073 ,0.274027 ,0.277225 ,0.284345 ,0.283710 ,0.283365 ,0.283367 ,0.286163 ,0.288472 ,0.288751 ,0.289256 ,0.290034 ,0.291031 ,0.288455 ,0.269891 ,0.199981 ,0.130150 ,0.127476 ,0.124950 ,0.122544 ,0.120275 ,0.118122 ,0.116089 ,0.114136 ,0.112292 ,0.110529 ,0.108873 ,0.107295 ,0.105787 ,0.104358 ,0.102985 ,0.101684 ,0.100442 ,0.099276 ,0.098167 ,0.097105 ,0.096074 ,0.095112 ,0.094205 ,0.093343 ,0.092513 ,0.091723 ,0.090978 ,0.090276 ,0.089609 ,0.088972 ,0.088379 ,0.087822 ,0.087297 ,0.086801 ,0.086339 ,0.085912 ,0.085511 ,0.085135 ,0.084798 ,0.084477 ,0.084196 ,0.083929 ,0.083712 ,0.083497 ,0.083322 ,0.083170 ,0.083040 ,0.082939 ,0.082867 ,0.082811 ,0.200259 }; TCanvas *c1 = new TCanvas("c1", "",71,231,600,600); gStyle->SetOptStat(0); gStyle->SetOptTitle(0); gStyle->SetHatchesLineWidth(1); gStyle->SetHatchesLineWidth(1); c1->SetFillColor(10); c1->SetBorderMode(0); c1->SetBorderSize(2); c1->SetTickx(); c1->SetTicky(); c1->SetLeftMargin(0.18); c1->SetRightMargin(0.07); c1->SetTopMargin(0.08); c1->SetBottomMargin(0.18); c1->SetFrameLineWidth(2); c1->SetFrameBorderMode(0); c1->SetFrameLineWidth(2); c1->SetFrameBorderMode(0); TH1F *hall = new TH1F("hall","",10,-5,95); TH1F *htpc = new TH1F("htpc","",90,0,90); TH1F *hset_etd = new TH1F("hset_etd","",90,0,90); TH1F *hvxd = new TH1F("hvxd","",90,0,90); TH1F *hsit_ftd = new TH1F("hsit_ftd","",90,0,90); // for(int i=0; i<90; i++)hall->SetBinContent(90-i,all[i]); for(int i=0; i<90; i++)hvxd->SetBinContent(90-i,vxd[i]); for(int i=0; i<90; i++)hset_etd->SetBinContent(90-i,set_etd[i]); for(int i=0; i<90; i++)hsit_ftd->SetBinContent(90-i,ftd_sit[i]); for(int i=0; i<90; i++)htpc->SetBinContent(90-i,tpc[i]); int fills = 3013; hallx = new TH1F("hallx","",10,0,90); hallx->GetXaxis()->SetTitle(""); hallx->GetYaxis()->SetTitle(""); hallx->SetMaximum(0.5); hallx->SetMinimum(0.); hallx->GetXaxis()->SetBinLabel(1,""); hallx->GetXaxis()->SetBinLabel(2,""); hallx->GetXaxis()->SetBinLabel(3,""); hallx->GetXaxis()->SetBinLabel(4,""); hallx->GetXaxis()->SetBinLabel(5,""); hallx->GetXaxis()->SetBinLabel(6,""); hallx->GetXaxis()->SetBinLabel(7,""); hallx->GetXaxis()->SetBinLabel(8,""); hallx->GetXaxis()->SetBinLabel(9,""); hallx->GetXaxis()->SetBinLabel(10,""); hallx->GetXaxis()->SetNdivisions(0); hset_etd->SetMaximum(0.5); hset_etd->SetLineColor(4); hset_etd->SetLineWidth(2); hset_etd->SetFillColor(4); hset_etd->SetFillStyle(fills); hall->GetXaxis()->SetTitle("#theta/degrees"); hall->GetYaxis()->SetTitle("X_{0}"); hall->SetMaximum(0.5); hall->SetMinimum(0.); hall->GetXaxis()->SetBinLabel(1,"90"); hall->GetXaxis()->SetBinLabel(2,"80"); hall->GetXaxis()->SetBinLabel(3,"70"); hall->GetXaxis()->SetBinLabel(4,"60"); hall->GetXaxis()->SetBinLabel(5,"50"); hall->GetXaxis()->SetBinLabel(6,"40"); hall->GetXaxis()->SetBinLabel(7,"30"); hall->GetXaxis()->SetBinLabel(8,"20"); hall->GetXaxis()->SetBinLabel(9,"10"); hall->GetXaxis()->SetBinLabel(10,"0"); hall->GetXaxis()->SetLabelSize(0.09); hall->GetXaxis()->SetTitleOffset(1.2); hall->GetXaxis()->SetNdivisions(209); htpc->SetLineColor(3); htpc->SetLineWidth(2); htpc->SetFillColor(3); htpc->SetFillStyle(fills); hsit_ftd->SetLineWidth(2); hsit_ftd->SetLineColor(2); hsit_ftd->SetFillColor(2); hsit_ftd->SetFillStyle(fills); hvxd->SetLineWidth(2); hvxd->SetLineColor(7); hvxd->SetFillColor(7); hvxd->SetFillStyle(fills); hallx->Draw(); hset_etd->Draw("SAME"); htpc->Draw("SAME"); hsit_ftd->Draw("SAME"); hvxd->Draw("SAME"); hall->Draw("AXIS SAME"); TLegend *leg = new TLegend(0.25,0.6,0.5,0.85,NULL,"brNDC"); leg->SetBorderSize(0); leg->SetLineColor(1); leg->SetLineStyle(1); leg->SetFillColor(10); leg->SetFillStyle(0); leg->SetTextSize(0.045); TLegendEntry *entry=leg->AddEntry("hset_etd","SET + ETD","F"); entry->SetLineStyle(1); entry->SetLineWidth(2); entry->SetFillStyle(fills); entry->SetLineColor(4); entry->SetFillColor(4); TLegendEntry *entry=leg->AddEntry("htpc","TPC","F"); entry->SetLineStyle(1); entry->SetLineWidth(2); entry->SetFillStyle(fills); entry->SetLineColor(3); entry->SetFillColor(3); TLegendEntry *entry=leg->AddEntry("hsit_etd","SIT + FTD","F"); entry->SetLineStyle(1); entry->SetLineWidth(2); entry->SetFillStyle(fills); entry->SetLineColor(2); entry->SetFillColor(2); TLegendEntry *entry=leg->AddEntry("hvxd","VTX","F"); entry->SetLineStyle(1); entry->SetLineWidth(2); entry->SetFillStyle(fills); entry->SetLineColor(7); entry->SetFillColor(7); leg->Draw(); c1->SaveAs("radLength.eps"); c1->SaveAs("radLength.gif"); }